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Abstract 

Stringent tests on top quark production and decay mechanisms are provided by the measure- 
ment of the top quark and W boson polarization. This paper presents a detailed study of these 
two measurements with the ATLAS detector, in the semileptonic (tt WWbb lyjijibb) 
and dileptonic {tt — > WWbb Ivlvbb) tt channels. It is based on leading-order Monte Carlo 
generators and on a fast simulation of the detector. A particular attention is paid to the system- 
atic uncertainties, which dominate the statistical errors after one LHC year at low luminosity 
(10 fb~^), and to the background estimate. Combining results from both channel studies, the 
longitudinal component of the W polarization (Fq) can be measured with a 2% accuracy and the 
right-handed component {Fr, which is zero in the Standard Model) with a 1% precision with 
10 fb~^ Even though the top quarks in ?f pairs are not polarized, a large asymmetry is expected 
within the Standard Model in the like-spin versus unlike-spin pair production. A 4% precision 
on this asymmetry measurement is possible with 10 fb^, after combining results from both 
channel studies. These promising results are converted in a sensitivity to new physics, such 
as tWb anomalous couplings, top decay to charged Higgs boson, or new 5-channels (heavy 
resonance, gravitons) in tt production. 
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1 Introduction 



Because of its high mass, intriguingly close to the electroweak symmetry breaking scale, the top 
quark raises interesting questions and its sector is an ideal place to look for new physics ElEI- 
Consequently, the search for non Standard Model interactions both in top quark production and 
decay is one of the main motivations for top quark physics. A consequence of the very high 
top mass is that this quark decays before it can form hadronic bound states 1,^. This unique 
feature among quarks allows direct top spin studies, since spin properties are not washed out 
by hadronization and since the typical top spin-flip time is much larger than the top lifetime. 
Therefore, top spin polarization f4l and correlation |5| are precisely predicted by the Standard 
Model (SM) and reflect fundamental interactions involved in the top quark production and de- 
cay. By testing only the top decay, the W boson polarization measurement complements top spin 
studies, helping to disentangle the origin of new physics, if observed. Namely, the t — > W^b 
decay mode is responsible for 99.9% of top quark decays in the SM. Therefore, the W polar- 
ization in the top decay is unambiguously predicted by the SM and its measurement provides a 
direct test of the tWb vertex understanding and more particularly of its V-A structure 16113. 

As a consequence, W polarization in top decay and top spin observables are sensitive 
probes of new physics in top production and decay. At the production level, a non-exhaustive list 
involves either anomalous couplings (HI El, which naturally arise in dynamical electroweak 
symmetry breaking models [2| such as technicolor 1 10| or topcolor 1 1 1 1, or new interactions, as 
for example a strong coupling of the top quark with a heavy spin resonance, such as a heavy 
(pseudo)scalar Higgs boson ifT^ as predicted e.g. by SUSY models ( gg ^ H ^ tt ), or the 
presence of extra dimensions lfT3l . At the decay level, deviations from the Standard Model can 
for example arise from tWb anomalous couplings, such as a V-i-A contribution in the vertex 
structure [il4il . or from a decay to charged Higgs boson lfT51 . 

Precise measurements of W and top polarization require a higher statistics than currently 
available from Tevatron data. The Large Hadron Collider (LHC) will be a top factory, pro- 
ducing more than 8 millions of tt events per year during its low luminosity running phase 
(10^-^ cm^^s^^), corresponding to an integrated luminosity of 10 fb^^. The production will 
occur through the gg tt (90%) and qq tt (10%) hard processes. Depending on the W 
decay modes, the tt events can be classified into three channels: the semileptonic channel 
{tt WWbb Zv jxjibb), the dileptonic channel {tt WWbb Ivlvbb) and the all hadronic 
channel {tt — > WWbb jijijsjAbb). The latter will be difficult to extract from the huge QCD 
background and has not been considered in this work. The electroweak single top production 
processes, which amount to approximately one third of the ti cross-section, can also be used to 
measure the polarization effects, but with a lower precision I.16J . They are not investigated in 
the following. 

The goal of this paper is to evaluate the precision to which W boson and top quark po- 
larization can be measured with the ATLAS detector, by combining results from dileptonic and 
semileptonic tt channels. It is organized as follows. Sectional discusses the W boson and top 
quark polarization in tt events, and gives the related physics observables. Section |3] presents 
the event simulation, reconstruction and selection, as well as a detailed background estimate. 
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Section |4l gives the expected ATLAS sensitivity to the W polarization, including a complete 
study of the systematic uncertainties. From these results, the sensitivity to the magnitude of 
tWb anomalous couplings that parametrize new physics is also extracted. Using the same se- 
lected events, section |5]presents the expected ATLAS sensitivity to the top polarization, and to 
the related physics beyond the SM. Sectionals dedicated to conclusions. 

2 W boson and top quark polarization in tt events 

This section presents the observables used to measure the polarization of the W boson (sec- 
tion |2]T]) and of the top quark (section l2!2b . 

2.1 W polarization observables 

The real in the t — > W'^b decay can be produced with a longitudinal, left-handed or right- 
handed helicity as sketched in Figure [T] The corresponding probabilities are Fq, Fi and Fr, 
respectively, whose SM expectations at tree level in the zero Z^-mass approximation are: 

' ^0 = ^^5^ = 0.703 + 0.002 x(M,- 175) 

< Ft = = 0.297 - 0.002 x (M^ - 175) (1) 

^ Fr = 0.000 

where Mf and Miy are the top and W masses in GeV. Left and right components are inverted for 
W bosons. By definition, we have the restriction F0+F1 + Fr = 1. Since massless particles 
must be left-handed in the SM, right-handed bosons do not exist in the zero Z?-mass ap- 
proximation, due to angular momentum conservation sketched in Figure [TJ Including QCD and 
electroweak radiative corrections, finite width corrections and non-zero b-quark mass induces 
small variations: Fq = 0.695, Fl = 0.304 and Fr = 0.001 for M, = 175 GeV O- Because the 
top quark is very heavy, Fq is large and the top decay is the only significant source of longitu- 
dinal W bosons ^ Deviations of Fq from its SM value would bring into question the validity of 
the Higgs mechanism of the spontaneous symmetry breaking, responsible for the longitudinal 
degree of freedom of the massive gauge bosons. Any deviation of Fr from zero could point to a 
non-SM V-i-A admixture to the standard left-handed weak current, as for example predicted by 
SU{2)l xSU{2)rxU{1) extensions of the SM [|18l . 

The best way to access particle spin information is to measure the angular distribution of 
its decay products, thereby called spin analyzers. As an example, illustrated in Figure |21 the 
charged lepton from the decay of longitudinally polarized tends to be emitted transversally 
to the direction, due to angular momentum conservation. Similarly, the charged lepton from 
a left-handed (right-handed) is preferentially emitted in the opposite (same) direction, 
leading to a softer (harder) p^ spectrum with respect to the leptons from longitudinal W^. The 

'QCD production, the only other source of real W bosons, produces nearly all W transversely polarized. 
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Figure 1: Sketches of angular momentum conservation in t — * W^b decay in the top rest frame. 
Simple (open) arrows denote particle direction of motion (spin). As a massless b-quark must be 
left-handed, the rightmost plot is forbidden in the SM at tree level. 

resulting angular lepton distributions are therefore very distinct for each W helicity state. 

As it is necessary to know the weak isospin of the W spin analyzer, the charged lepton 
is the best choice since w-like jets can not be distinguished experimentally from J-like jets. 
Consequently, the W polarization is better measured in dileptonic and semileptonic tt channels 
through the distribution of the m angle between the charged lepton direction in the W rest frame 
and the W direction in the top quark rest frame. The ^ angular distribution is given by the 
following expression |l6||: 



1 dN 



Fo 



/sin^Y 

Vvt) 



1-cos^ 



+ Fr 



1 +cos^ 



(2) 



Its SM expectation is shown in Figure |3l It reflects the superposition of the three terms of 
Equation Q, corresponding to the longitudinal (sin^)^, the left-handed (1 — cos^)^ and the 
right-handed (1 + cos^)^ helicity states. Each term is weighted by the fraction Fq, Fl or Fr 
given in Equation (HJ). 
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Figure 2: Sketches of the different polarization modes in t W^b decay and resulting 
lepton directions. Simple (open) arrows denote particle direction of motion (spin). For W , left 
and right-handed components are inverted. 




Figure 3: Angular distribution of Equation in the SM. The predicted contributions from 
longitudinal (0) and left-handed (L) helicity states are shown separately with dashed lines. The 
right-handed contribution is null in the SM. The sum (0+L+R) is depicted with a full line. 
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Since the W and top rest frames are used in the ^ angle measurement, it requires a com- 
plete event topology reconstruction. This is rather easy in the semileptonic tt channel, with only 
one neutrino in the final state and a high signal over background ratio (see section 17.4. II) . In 
the dileptonic channel^, the event reconstruction is more challenging (see section lT.4.21 ). There- 
fore, the ^ angle is reconstructed in terms of the invariant mass of the lepton and the Z?-quark, 

cos^ , '\ -1 (3) 

which is valid in the zero b-mass limit, and where M/ and Mw are set to 175 GeV and 80.41 GeV, 
respectively. In this approach, the dependence on the Z7-jet energy scale and on the top mass un- 
certainty is high. On the contrary, these two dependences cancel at first order when measuring 
directly cos^. A study performed in the semileptonic channel shows a two times lower sys- 
tematics on Fq, Fl and Fr using the cos^ observable compared to that obtained with M^^. 

2.2 Top polarization observables 

Similarly as for the W , the top polarization can be analyzed trough the angular distribution of 
its daughters. In this case, the spin analyzer, denoted by z, can be either a direct daughter (W, b) 
ox aW decay product (/, v, ji or J2). The relevant angular distribution is lfT4ll : 



^7^ = i(l+^«'^«^^') (4) 
N a cos 0/ 2 

where S is the modulus of the top polarization and 0/ is the angle between the direction of 
particle i in the top quark rest frame and the direction of the top polarization, a, is the spin 
analyzing power of this particle. It is the degree to which its direction is correlated to the spin 
of the parent top quark. It has been computed at the next-to-leading order (NLO) since long 
for the lepton (/) ll23ll and more recently for the b quark, the W boson and the quarks from the 
W decay Il24ll . The theoretical values are given both at LO and NLO in Table [T] for a spin up 
top quark (signs are reversed for a spin down or for an anti-quark). Even if W and b are direct 
daughters of the top, their analyzing power is low due to the intrinsic polarization of the W 
which interferes destructively with the top one. Consequently, at LO, the difference Fq — Fi 
provides a measurement of ttiy |25|. Charged leptons and down-type quarks, which are almost 
100% correlated with the top spin direction, are optimal spin analyzers. But on the contrary 
to leptons, d and 5-jets cannot be distinguished experimentally from u and c-jets. Therefore, 
the analyzing power of light jets is the average value ttjet ~ (1 — 0.31) /2 = 0.35. This can be 
improved by choosing the least energetic jet {le j) in the top rest frame, which is of d type in 
61% of the case, resulting to a, ~ 0.5 ll24l . 

Equation Q can be directly used for top quarks produced lonely via the weak interaction, 
which are polarized (5 ~ 1). This is not the case of the top quarks produced in tt pairs, which are 
not polarized lE^ . However, the top and the anti-top spins are correlated, which can be easily 
understood. Close to the tt production threshold, the tt system produced by qq annihilation is 

^Tliis is also the case at tlie Tevatron 1 191 1201 12 II 1221 and for single top analysis at LHC 1161 . 
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Particle 


Z?-jet 


w+ 


/+ 


ji=J-jet, 5-jet 


j2=M-jet, c-jet 


lej 


a, (LO) 


-0.41 


0.41 


1 


1 


-0.31 


0.51 


a, (NLO) 


-0.39 


0.39 


0.998 


0.93 


-0.31 


0.47 



Table 1: SM spin analyzing power at LO and NLO of top quark daughters: b-jet, W^, and W 
decay products : lepton (Z+j, ji, j2 or the least energetic non b-jet in the top rest frame, called 
le j lt24J . The top quark is spin up. Signs are reversed for a spin down or for an anti-quark. 



in a ^Si state, while it is in a ^5o state with the gg fusion process. Therefore, in the first case, 
the top quarks tend to have their spins aligned while in the second case their spins tend to be 
opposite to each other. Away from threshold, this simple picture is modified due to the presence 
of angular momentum. In the absence of polarization, a direct measurement of the correlation 
at the level of the top quarks is obtained from observables of the form: 

(a-S,)(b-S,-), (S,-S,-) (5) 

where St, St are the spin operators of the top and anti-top, and a, b are arbitrary directions 
(|a| = |b| = 1). A more familiar representation of the observables shown in Equation ^ can be 
obtained from the relation: 

A^4((a.s,)(s.s,)>= "f'^;+"f^^;-"w;-"f^!^; (« 

<5(<,<,) + <5(<j<j) + <5(<,<j) + 0(<j<,) 

where o(?| /|f| /|) denotes the cross section for the production of a top quark pair with spins up 
or down with respect to a quantization axis defined by a in case of the top quark and b in case 
of the anti-top quark. Note that rewriting the observable (S/ • Sf) using 

Ao = {St ■ St) = Y^St^Sn = £(e(') ■ S,) (e(') ■ S,-) (7) 

i i 

where e['^ = ht^, the observable {St ■ S-f) can also be cast into the form shown in Equation Q. 
Computation in the Standard Model gives favorably like-spin pairs (A > 0) when using the 'he- 
licity' basis for the spin basis at LHC |15|. In this basis-^, the top (anti-top) spin quantization 
axis corresponds to the top (anti-top) direction of flight in the tt center of mass system, and the 
notation i is replaced by L (Left) and R (Right). Figured shows the invariant mass distribu- 
tion of the tt system with like and unlike helicities for the two possible production mechanisms. 
As already explained, gg and qq processes contribute to the asymmetry with opposite signs 
(A > for gg and A < for qq). The theoretical Standard Model value integrated over the 
whole LHC spectrum at LO is A = 0.319 and Az)=-0.217. At NLO these values become i^: 

A = 0.326t^Z{^)X^l\iPDF) 

(8) 

Ao = -0.23lt^Z{,)t^ZiPI^F) 

Systematic uncertainties come from factorization and renormalization scales (p = /jf = ^r) and 
from Parton Distribution Function (PDF). As A and A/j are defined as ratios between two cross 



^Another basis was recently found to be more optimal, but more complicated to reconstruct l27l . 



6 



sections, PDF, and a5 dependences cancel to a large extent. Moreover, NLO QCD corrections 
are small and thus theoretical uncertainties are well under control. 

At LHC, it is possible to increase the asymmetry by applying an upper cut on the tt invari- 
ant mass Mff. As shown in Figured the asymmetry is maximal at low invariant masses for the 
gg contribution, which is by far dominant at the LHC, and is equal to around M^f = 900 GeV. 
Therefore, selecting low energetic top quarks with Mtf < 550 GeV rejects only 30% of the 
events while A and Ao are enhanced by about 30% at LO: 




400 600 800 1000 1200 1400 1600 , 

M„ (GeV/c^) 



Figure 5: Up, double differential angu- 
Figure 4: Invariant mass distribution of the tt sys- lar distribution of Equation (|70l). Down, 
tern with like (LL+RR) and unlike (LR+RL) helic- opening angle distribution of Equa- 
ities for the two possible production mechanisms tion (|77I- Both are shown at parton 
(gg and qq) level with M,t < 550 GeV. 
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Similarly to Equation Q for top polarization, angular distributions can be used to probe 
the tt spin correlation, as: 



the double differential angular distribution of top and anti-top quark decay products ll30ll : 
1 (fN 1 



-(1 +5iCOS0i +fi2COS02 -CCOS01COS02) (10) 



A^<icos0i<icos02 4 

where 0i (02) is the angle between the direction of the t (t) spin analyzer in the t (f) rest 
frame and the t (J) direction in the tt center of mass system. These complicated angles are 
the direct consequence of the choice of the 'helicity' basis for measuring the asymmetry. 
As top and anti-top quarks are not polarized in this basis, B\ = B2 = 0. Figure |5l (top 
panel) illustrates this double angular distribution for the SM in the semileptonic channel. 

the opening angle distribution lE91l : 

1 ^'l-Dcos4>) (11) 



Ndcos^ 2 

where 4> is the angle between the direction of flight of the two spin analyzers, defined 
in the t and t rest frames respectively. Figure |5l (down panel) shows this opening angle 
distribution for the SM in the semileptonic channel. 

In Equations (flOb and (fTTT) . C and D are the spin correlation observables. Before any 
phase-space cut, they can be easily measured using the following unbiased estimators [i31il : 

C — —9 < cos 01 cos 02 > 

(12) 

D = -3 < cos4> > 

The production asymmetries A and Ad are then directly deduced by simply unfolding the decay 
contribution through the spin analyzing power of the daugther particles tti and a2: 

A- c 



Ad 



aia2\ 

(13) 

D 

' I«ia2| 



Since leptons are the most powerful spin analyzers (a = 1), the dileptonic case is a priori the 
most promising. On the contrary, the all hadronic case is the most unfavorable, with a low 
spin analyzing power and a huge background. Several choices exist in the semileptonic channel 
for the spin analyzers on the hadronic side (W, b and lej) and experimentally the lej is the 
best choice ll32l l. In any case the semileptonic channel is more challenging compared to the 
dileptonic one because its spin analyzers are less powerful than the leptons. But on the contrary 
to the dileptonic channel, the number of events is 6 times larger and the event reconstruction is 
much easier (only one neutrino in the final state). 
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3 Event simulation, selection and reconstruction 



This section describes the software tools used to generate and simulate signal and background 
events using a modeling of the ATLAS detector. Then the event selection and reconstruction is 
explained in both tt semileptonic and dileptonic channels. 

3.1 Signal and background definition 

General figures of ?f pairs decay are: Br{t Wb) ~ 1, Br{W Zv/) ~ 1/3 with I = x in 
equal probabilities and Br{W — > qiqi) ~ 2/3 with qi^qi) = u{d),c{s) in equal probabilities. 
With a NLO cross-section around 850 pb ll33ll . 3.8 (0.9) millions of tt semileptonic (dileptonic) 
events will be produced with an integrated luminosity of 10 fb ^ corresponding to one LHC 
year at low luminosity. Among them, 2.5 (0.4) millions are signal events, defined as: 



with Z = The 1.3 (0.4) million events with at least one Z = x in the semileptonic (dileptonic) 
channel are considered as background. Non-?f background is composed of QCD background, 
from which mainly bb production is relevant for our study, and of electroweak backgrounds, 
which are V7+jets, Z{-^ ZZ)+jets, Wbb, 2 vector bosons (ZZ, ZW , WW) and single top produc- 
tion. 

3.2 Event generation 

The Monte Carlo leading-order generator TopReX 4.05 ['341 is used for the tt event generation. 
It includes the Standard Model LO tt spin correlation^. A top mass of 175 GeV is assumed 
and the structure function CTEQ5L is used. The 2^-scale (pr(?)^ -fM^) used for a5 is the 
same as for the structure function. The proportion of gg and qq processes, which directly im- 
pacts the spin correlation (see Figure |4|), is 86%/ 14%. Partons are fragmented and hadronized 
using PYTHIA 6.2 |37|, including initial and final state radiations, as well as multiple interac- 
tions, in agreement with CDF data extrapolated to LHC |38 1. The Z^-fragmentation is performed 
using the Peterson parametrization with = —0.006. TAUOLA and PHOTOS [39| are used to 
process the x-decay and radiative corrections. All results correspond to one LHC year at low 
luminosity. For the systematics study, samples corresponding to three times (ten times) more 
statistics are generated for each source of uncertainties in the semileptonic (dileptonic) channel. 

For what concerns the non-?f background generation in the semileptonic channel, PYTHIA 
is used, except for W+A jets and Wbb which are treated with AlpGen [40| and AcerMC [ 4T]| 
generators, respectively. About 3 • 10^^ W{-^ ly)+A jets weighted events are generated with 
cuts on the four extra light jets: pr >10 GeV, |ri| < 2.5 and Ai?(jet-jet) >0.4^. Despite this 
huge effort, it only represents 1/63^^^ of the statistics for one LHC year (380000 events). For bb 

^NLO spin correlation simulations are expected to be included in the future in MCatNLO 1351 generator. 
^A/; = v/A(|)2+Ar|2 




(14) 
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background, given the very high cross-section (~500 ^ub) of the process, the cut >120 GeV 
is applied at the parton level. 750 million of events have been generated, corresponding to 1/8'^'' 
of the statistics for one year. Anyway, as this QCD background is very difficult to estimate, it 
should be extracted from the data. Except for W+4 jets and bb, the statistics corresponding to 
one LHC year at low luminosity is simulated for each background without any cut at the parton 
level. In the dileptonic case, all backgrounds are simulated using PYTHIA 6.2. 



3.3 Detector modeling 

A simplified modeUng of the ATLAS detector, ATLFAST 2.6.0 ||42|, is used. This essentially 
accounts for resolution smearing of objects accepted within the detector geometry, according to 
the expected performances |43 1. Only settings of particular importance are recalled here: 

• Isolation criteria only for lepton (electron and muon) consists in: /) asking Ej < 10 GeV 
in a pointing cone of 0.2 around the lepton and ii) requiring the nearest calorimeter cluster 
at Ai? > 0.4. 

• Jets are reconstructed with a cone algorithm, with a size AR = 0.4. They are calibrated to 
obtain a correct jet energy scale 1.43 J . 

• A 60% Z^-tagging efficiency is assumed, as well as a c-jet rejection of 10. For the other 
jets, the rejection is 100. This is a rather pessimistic assumption compared to the latest 
simulation results 1 44 1 . 

No trigger inefficiencies and no detailed acceptance (as crack between barrel and endcap) are 
included in this analysis. A 90% lepton and a 95% jet reconstruction efficiency are assumed. 

3.4 Event selection and reconstruction 

The heavy top mass makes the event topology at LHC outstanding: t and f are preferentially 
produced in the central region (|ri| < 2.5) and back to back in the transverse plane. Therefore, t 
and t are naturally well separated, which will ease the event reconstruction. 

3.4.1 Semileptonic tt events 

Semileptonic signal events are characterized by one (and only one) isolated lepton, at least 4 jets 
of which 2 are b-jets, and missing energy. They are selected by requiring an isolated lepton from 
first and second level trigger. The offline kinematic cut on the lepton is directly given by the 
trigger threshold which is set to pr > 20 GeV at 10-^^cm^^s^^ Moreover, at least four jets 
with pt > 30 GeV are required, among which at least two are Z7-tagged. This pr cut is a good 
compromise between a low combinatorial background and a good statistics. The 20 GeV cut on 
the missing transverse energy (p^"'') is standard for events with one neutrino and rejects almost 
no signal events. All these kinematic cuts are summarized in Table|2l Their resulting efficiency 
on signal events is 6.5%. The average pr of the lepton and least energetic jet in the top rest 
frame {le j) are around 50 GeV. 
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After kinematic cuts, the full event topology is reconstructed. On the top hadronic decay 
side, the two non Z^-jets with Mjj closest and within 20 GeV to the known M\y [45 1 are selected. 
Then, the b-jet with Mjjt, closest to the known M, P31 is chosen to reconstruct the hadronic top 
and removed from the Z7-jet list. On the leptonic side, the missing transverse momentum is used 
to evaluate the neutrino p^. Its longitudinal component, p^, is determined by constraining M/v 
to Miv: 



+ ^{pW + iPlY - {P[+Plf - iPy+P'y)' - {P[+Ply 



M 



w 



(15) 



For 25% of the events, there is no solution since p™"'^ overestimates p^. In these cases, p^ is 
decreased step by step by 1% until a solution is reached f46|. Then, for the 5% of events with 
more than two Z^-jets, the one closest to the lepton in A/? is chosen to reconstruct the leptonic top. 

Figurel^(a-c) shows reconstructed W and top masses. Results are comparable with those 
of the top mass study PTll . Quality cuts are then applied on top and anti-top reconstructed 
masses (|Mf°'^— M/| < 35 GeV and |m'^'^— M,| < 35 GeV) to reject badly reconstructed events. 
At this stage, 3.3% of the signal events are kept, corresponding to 85000 signal events for one 
LHC year at low luminosity. TableEllists all the selection cuts. 



Selection type 


Variables 


Cuts 




=1 isolated lepton 


pr > 20 GeV, |ri| < 2.5 


Kinematic and 


> 4 jets 


pr > 30 GeV, ri < 2.5 


acceptance 


Z?-tagged jets 


> 2 




Missing energy (v) 


p^'"- > 20 GeV 




|M'-^-Miv| 


< 20 GeV 


Reconstruction 


|Mf"^-Mf| 


< 35 GeV 


quality 




< 35 GeV 



Table 2: Selection cuts in the semileptonic tt channel. 



After selection criteria are applied, the background is composed for more than 80% of 
tt ^ x + X events, as shown in Tabled The amplitude and shape of this ?f background should 
be easily under control. The remaining non-?f background is dominated by /v)+4 jets, bb 
and single top events. In the first two cases, since only a few tens of events subsist after all cuts, 
Poisson statistics is used to give an estimate of the expected number of events. Given its very 
low contribution to the overall background, the non-tt background will be neglected in the rest 
of the analysis. In total, about 7000 background events are expected for one LHC year at low 
luminosity, giving a signal over background ratio of 12. 

For the top spin study, to enhance the correlation, a further cut, M^f < 550 GeV, is applied 
on the tt reconstructed mass (section IT^ . whose distribution is shown in Figure |^(d). The 
total efficiency becomes 2.3%, corresponding to 60000 signal events for one LHC year at low 
luminosity. 
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Figure 6: (a-c): Reconstructed masses of W from hadronic side, top from hadronic and leptonic 
sides in the semileptonic channel. Lines correspond to Gaussian+third order polynomial fits, 
from which Mean and o are extracted, (d): Normalized reconstructed mass distributions of the 
tt system are shown both in the semileptonic (full) and dileptonic (dashed) channels. 





Initial number of 


Number of 




Events (xlO^) 


selected events 


Signal (tt semileptonic) 


2.5 


85000 


tt background 






tt-^x+X 


1.3 


6200 


tt all had 


3.7 


70 


Non-tt background 






/v)+4jets 


24 


[400,1000] 


bb (v^>120 GeV) 


6000 


200 


Z(-^ Z/)+jets 


49 


12 


zz,ww,zw 


1.1 


5 


l\)bb 


0.7 


3 


single top 


1.0 


350 



Table 3: Number of events in the semileptonic tt channel (signal and background) for one LHC 
year at low luminosity, 10 fb^^, before and after selection cuts. 
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3.4.2 Dileptonic tt events 

Dileptonic events are characterized by two (and only two) opposite charged isolated leptons, at 
least two jets of which two are Z?-jets, and missing energy. They are selected by requiring two 
leptons from first and second level trigger. The offline p^ cut on opposite sign leptons is con- 
servatively set to 20 GeV, well above the trigger thresholds which are lower or equal to 15 GeV 
for two leptons. Moreover, two Z?-tagged jets with p^ > 20 GeV are required. The 40 GeV 
cut on p^"'' is standard for events with two neutrinos. Tabled shows the selection cuts in the 
dileptonic channel. Their total efficiency on signal events is 6.5%. 



Selection type 


Variables 


Cuts 


Kinematic and 
acceptance 


=2 isolated leptons 
> 2 jets 
^-tagged jets 
Missing energy (v) 


pr > 20 GeV, |ri| < 2.5 
pr > 20 GeV, \r\\ < 2.5 
= 2 

pf'' > 40 GeV 



Table 4: Selection cuts in the dileptonic tt channel. 



After kinematic cuts, the event topology is reconstructed using the algorithm developed 
in EHJ. The aim of the reconstruction is to obtain the unknown momenta of neutrino and anti- 
neutrino and the association between the two Z?-jets and the b and b quarks. To solve the set of 
six non-linear equations coming from the momenta and energy conservation, the known M^ and 
Miy are assumed. The set of equations can have up to four solutions for each combination of 
the association b-jets to b and b quarks. The choice of the solution is based on the computation 
of weights from known distribution of transversal momenta of t, t and v, v. The reconstruction 
efficiency of this algorithm is 80% with the correct solution found in 65% of the cases. Most of 
the dilution comes from the wrong b assignation. After cuts and reconstruction, 5.3% of the sig- 
nal events are kept, corresponding to 21000 signal events for one LHC year at low luminosity. 
The background is then composed for 90% oitt ^ x + X events, as shown in Table |5l In total, 
less than 4000 background events are expected for one LHC year at low luminosity, giving a 
signal over background ratio of 6, two times lower than in the semileptonic channel. 

As in the semileptonic case, a cut on the tt reconstructed mass, M^^ < 550 GeV, whose 
distribution is shown in Figure|6l(d), dashed lines, is applied to enhance the spin correlation. In 
this case, the total efficiency becomes 3.5%, corresponding to 15000 signal events for one LHC 
year at low luminosity. 
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Number of 


Number of 




Events (x 10^) 


selected events 


^inrnol (ft f\^\f^r^\f^\T\^n\ 




91000 
z i www 


tt background 






tt — > %-\- L 


U.J 


J /uu 


tt^l + jet 


3.8 


40 


Non-tt background 






bb ipT >20 GeV) 


30000 


<200 


Z+jets, W+jets, ZZ, WW,ZW,Wbb 


4500 


<100 


single top 


1.0 


7 



Table 5: Number of events in the dileptonic tt channel (signal and background) for one LHC 
year at low luminosity, 10 fb^^, before and after selection cuts. 

4 Sensitivity to W boson polarization in tt events 

In this section the method to extract the W polarization observables is explained ("section 1431) . 
the complete systematics study is presented (section |421) and results combining dileptonic and 
semileptonic channels are given (section l431) . Finally, using these results, the sensitivity to tWb 
anomalous couplings is discussed (14.41) . 

4.1 Measurement method 

As charged leptons from longitudinal W have a harder p^ spectrum than those from left- 
handed W , they are more likely to pass the trigger threshold and offline selection requirements. 
More generally, the reconstructed cos^ angular distribution is distorted by all the acceptance 
and reconstruction effects compared to the parton level one, as shown in Figure IT] This leads 
to a bias in the measurement, which is more pronounced as the cut on the lepton pj- increases. 
To correct for it, a weight is applied on an event by event basis, allowing to recover, as much as 
possible, the original shape. The weighting function is obtained from the ratio between the two 
normalized distributions of Figure |7] (i.e. after selection cuts and at parton level) for semilep- 
tonic and dileptonic events. In both channels, this ratio, shown in Figure [HI is fitted by a third 
order polynomial function to extract a smooth correction. The fit is restricted to the region 
—0.9 < cos^ < 0.9, which is the most extended region where the correction is varying slowly. 
The correction functions, computed on an independent data sample, are then applied event by 
event on the analysis samples. 

The final distribution after event selection and correction is shown in Figure|9lfor semilep- 
tonic and dileptonic events, corresponding to one year at low luminosity, 10 fb ^ The W po- 
larization is extracted from a fit in the restricted region with the Equation dJJ) function and the 
constraint Fq-\- Fi-\- Fr = 1. The results for Fq, Fi and Fr are compatible with their SM ex- 
pectations. The statistical errors are 0.005 for Fq and 0.003 for Fi and Fr in the semileptonic 
case and increase to 0.010, 0.007 and 0.005 in the dileptonic case. The correlation between the 
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COS COS * 



Figure 7: Normalized cos^ distribution after selection cuts (histogram) for semileptonic (left) 
and dileptonic (right) tt events. For comparison, the SM parton level distribution (Figure^Sj) is 
superimposed (full line). 




cos cos 



Figure 8: Ratio between the two normalized distributions o/cos^ shown in Figure^ i.e. after 
selection cuts and at parton level, for semileptonic (left) and dileptonic (right) tt events. The 
full line is the result of a third order polynomial fit in the range —0.9 < cos^ < 0.9. 
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Figure 9: Normalized cos^ distribution after selection cuts and correction for semileptonic 
(left) and dileptonic (right) tt events after one year at low luminosity, 10 fb^^. The full line is the 
result of a fit with function of Equation in the range —0.9 < cos ^ < 0.9, with Fq + Fi + Fk = 
1. The dashed line represents the continuation of the function outside the fit region. 

parameters for the couples (Fo,Fi), (Fo,F^) and {Fi,Fr) are -0.9, -0.8 and 0.4. 

The systematic uncertainty induced by the weighting method has been estimated by vary- 
ing the number of bins (from 40 to 25), the fit limits (from [-0.9;0.9] to [-0.8;0.8]) and the 
polynomial order (from P3 to P5 and P7). All variations are below the statistical error. There- 
fore, the total uncertainty of the method is estimated to be smaller than the statistical error. 

The corrections functions of Figure [Hare extracted with a Standard Model scenario, as- 
suming a pure V-A top decay vertex. In case of deviation from the SM, the kinematic distribu- 
tions, such as lepton p^ or angles can be affected. This is for example the case if a V-i-A com- 
ponent is present. In such a scenario, the fraction of longitudinal W bosons will be unchanged, 
but a right component Fr, whose lepton spectrum is harder (see section lOT) will appear. As a 
consequence, the correction function will be changed. This is illustrated in Figure [TOlfor differ- 
ent Fr input values. In each case, the statistics of one LHC year of semileptonic events has been 
generated with AlpGen. Applying the SM correction function to these samples will therefore 
not correct completely for the bias induced by the selection cuts. Figure [TT] shows the fraction 
Fr extracted from the fit^ as a function of the Fr input value after applying the SM correction 
function (open circles). The measurement is clearly biased. To overcome this problem we pro- 
ceed iteratively. The SM correction function is first used. Then, in case of deviation of Fr from 
zero, a new correction function is calculated with this new Fr component, and applied. The 
process converges after a few iterations, as shown in Figure [TT1 

^In this case, Fq is fixed to its SM value and the only fitted parameter is Fr. 
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Figure 10: Correction functions for different 
Fr input values in the semileptonic channel. 
The case Fr = Q corresponds to the SM, and 
the case Fr = 0.297 corresponds to a pure 
V+A model. The SM function (pure V-A) is 
the one of Figure^ 
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Figure 11: Fr extracted from the fit as a func- 
tion of Fr input values. The open circles cor- 
respond to the measurement using the SM cor- 
rection. Full triangles, stars and squares cor- 
respond to the first, second and third iteration, 
respectively ( see text). The dotted line is y = x. 



4.2 Systematic uncertainties 

This section presents a detailed study of the systematic uncertainties related to the W polariza- 
tion measurement. 



4.2.1 Systematic uncertainties at generation level 

In this subsection, five main sources of systematic uncertainties are considered: 

• 2-scale: The uncertainty related to the Q-scale at parton generation is estimated by com- 
paring samples generated with TopReX and AlpGen using default 2^-scale: prit)^ + Mj 
and M^. In both cases, the same hadronization scheme (PYTHIA) is used. 

• Structure function : The impact of the structure function is estimated by the maximum 
difference between the measurements obtained with the standard parton density function 
(CTEQ5L) and three other ones, CTEQ6L |49|, MRST2002 [50| and GRV98 (Hll- It is 
conservatively estimated to be independent of the above Q-scale uncertainty. 

• ISR, FSR: The presence of initial state radiations (ISR) from incoming partons and espe- 
cially final state radiations (FSR) can affect the ^ angle reconstruction since it impacts 
the top quark reconstruction. To estimate the effect due to ISR, the difference between the 
measurements obtained with ISR switched on (usual data set) and off is computed. The 
same approach is used for FSR. The level of knowledge of ISR and FSR is around 10%. 
Therefore, as a more conservative estimate, the systematics uncertainties have been taken 
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to be 20% of the corresponding differences. It should be noted that more sophisticated 
methods exist to make this evaluation [52|. 

• Z7-fragmentation: The Z?-quark fragmentation is performed according to the Peterson para- 
metrization^, with one free parameter £/,. The default value is set to = —0.006. It has 
been changed to a more recent LEP value (£^ = —0.0035 [154] ). and the differences on the 
results are taken as systematic uncertainties. 

• Hadronization scheme : The angular distributions of jets and leptons may be influenced 
by the hadronization scheme. Generating partons with AcerMC then processing the 
hadronization with PYTHIA or HERWIG ll55l leads to different W polarization mea- 
surements. The related systematic uncertainty is extracted from this difference. 

4.2.2 Systematic uncertainties at reconstruction level 

In this subsection, three main sources of systematic uncertainties are considered: 

• Z?-tagging: The impact of the Z?-tagging efficiency is studied by increasing it from 50% 
to 70% by steps of 5%, according to a parametrization coming from full simulation iH^ . 
Increasing the Z?-tagging efficiency degrades the c-jets and light jets rejection factors. As 
an example, going from 55% to 60% decreases them respectively by 30% and 80%. Fig- 
ure [T2I (a) shows Fi, Fq and Fr measurements as a function of the Z?-tagging efficiency. 
Small and smooth dependences are observed. The related error is computed with a real- 
istic ±5% uncertainty on the Z?-tagging efficiency. 

• Z7-jet miscalibration: The impact of the knowledge of the absolute b-jet energy scale is 
estimated by miscalibrating the reconstructed b-jet energy. Results are shown in Fig- 
ure [T21(b), for a miscalibration between ±5% by steps of 2%. The behaviors can be 
easily understood: a positive miscalibration overestimates the invariant mass of the lep- 
ton and the Z?-quark, M/^ and therefore cos^. Equation ©. This bias the polarization 
toward higher values {Fi decreases and Fr increases). The corresponding systematic is 
computed with a realistic ±3% uncertainty on the Z7-jet energy scale. 

• Input top mass: The SM W polarization has a small dependence on the top mass, with 
an increase (decrease) of Fq (Fi) by 0.002 per GeV, Equation ([T]). Moreover, a 175 GeV 
top mass value is assumed in the event reconstruction, which can impact the cos^ mea- 
surement if the real top mass is different. Therefore, different samples of events were 
generated with a top mass between 170 and 180 GeV by steps of 2.5 GeV. Results are 
shown in Figure [T21(c). As for the case of a positive b-jet miscalibration, a high top mass 
increases Mit, and therefore bias the polarization toward higher values. The related sys- 
tematics are -0.008 (-0.010) per GeV on Fl, 0.006 (0.003) on Fq and 0.002 (0.007) on 
Fr in the semileptonic (dileptonic) channel. The final uncertainty is computed assuming 
AMtop=2 GeV, which should be reached at Tevatron run II ll5Sll . 



^It has to be noticed that recent measurements can not be well fitted with the Peterson parametrization II53I . 
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Figure 12: Measured Fi (left), Fq (middle) and F^ (right) in the semileptonic (black circles) 
and dileptonic (open circles) tt channels as a function of different parameters, see text for more 
details. Linear fits are superimposed in each case, and the corresponding slope is indicated. 
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4.2.3 Other sources of systematic uncertainties 

• Background: As shown in section l34l the only sizable background comes from tt^x+X 
events, and is well under control. A large variation of this background level by ±10% has 
a negligeable impact on the results when considering the semileptonic channel because 
of the high signal over background ratio. In the dileptonic channel, this results in an 
uncertainty of 0.004 on Fl, 0.003 on Fq and 0.001 on Fr. 

• Pile-up: Pile-up events may influence the reconstruction and therefore impact on W po- 
larization measurement. They are generated with PYTHIA setting MSTP(131)=1 and 
using the default process definition MSTP(132)=4. Different samples were generated 
adding 2.3 or 4.6 pile-up events according to the Poisson law. These numbers are ex- 
pected for a luminosity of 10-^^cm^^s^^ and 2x 10^-^cm^^s^^, respectively. Results are 
shown in Figure [121 (d). It can be noticed that these results are conservative, as no jet 
recalibration is applied. 

4.2.4 Systematics summary 

All systematic uncertainties are listed in Table |6l and illustrated in Figures [13] and [141 for the 
semileptonic and dileptonic channels, respectively. Generation and reconstruction sources con- 
tribute roughly in the same proportion to the total error. The dominant generation contributions 
come from the FSR knowledge, the hadronization scheme and the Q-scale. while the recon- 
struction systematics are dominated by the b-jet miscalibration and the top mass uncertainty. 



Source of uncertainty 


Semileptonic channel 


Dileptonic channel 




Fl 


Fo 


Fr 


Fl 


Fo 


Fr 


Generation 














2- scale 


0.000 


0.001 


0.001 


0.011 


0.010 


0.002 


Structure function 


0.003 


0.003 


0.004 


0.002 


0.004 


0.003 


ISR 


0.001 


0.002 


0.001 


0.000 


0.001 


0.001 


FSR 


0.009 


0.007 


0.002 


0.016 


0.008 


0.008 


Z^-fragmentation 


0.001 


0.002 


0.001 


0.002 


0.003 


0.002 


Hadronization scheme 


0.010 


0.016 


0.006 


0.002 


0.003 


0.002 


Reconstruction 














Z7-tagging (5%) 


0.006 


0.006 


0.000 


0.003 


0.003 


0.001 


b-jet miscalibration (3%) 


0.011 


0.005 


0.005 


0.017 


0.001 


0.017 


Input top mass (2 GeV) 


0.015 


0.011 


0.004 


0.021 


0.007 


0.014 


Others 














S/B scale (10%) 


0.000 


0.000 


0.000 


0.004 


0.003 


0.001 


Pile-up (2.3 events) 


0.005 


0.002 


0.006 


0.001 


0.002 


0.003 


TOTAL 


0.024 


0.023 


0.012 


0.034 


0.016 


0.024 



Table 6: Summary of systematics on Fl, Fq and Fr in the semileptonic and dileptonic tt channels. 
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Figure 13: Systematic uncertainties on Ft, Fq and Fr in the semileptonic tt channel. 
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Figure 14: Systematic uncertainties on Fl, Fq and Fr in the dileptonic tt channel. 
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4.3 Results 



Table III presents the expected Standard Model results for the W polarization measurement in 
tt semileptonic and dileptonic channels after one LHC year at low luminosity (lO-'^cm^^s^^ 
10 fb^^). The sensitivity is driven by the systematic uncertainties, which largely dominates the 
statistical ones. Fi and Fr parameters are more precisely measured in the semileptonic channel, 
while the accuracy on Fq is slightly better in the dileptonic one. Combining the results of both 
channel studies, assuming a pessimistic 100% correlation of systematic errors, lead to the results 
shown in Table|7l rightmost column. The only improvement of the combination concerns Fq on 
which the absolute error is estimated to be 0.016. It is worth to notice that Fr, which is expected 
to be zero in the SM, is the most precisely measured with an accuracy of 0.012. 





Semileptonic (±stat±syst) 


Dileptonic (±stat±syst) 


Semilep+Dilep 


Fl 
Fo 
Fr 


0.299 ± 0.003 ± 0.024 
0.699 ± 0.005 ± 0.023 
0.002 ± 0.003 ± 0.012 


0.314 ±0.007 ±0.034 
0.696 ±0.010 ±0.016 
-0.010 ± 0.005 ± 0.024 


0.303 ± 0.003 ± 0.024 
0.697 ± 0.004 ± 0.015 
0.000 ± 0.003 ± 0.012 



Table 7: Standard Model results for W polarization components after one LHC year of data 
taking (lO^^cmT^s^^, 10 fb^^) in semileptonic and dileptonic tt channels. A combination of 
both results is presented in the last column. 

This result is 3 times better than the statistical error foreseen with single top events (~ 0.03 
on Fr with 10 fb~^) lfT6ll . It is also roughly 3 to 5 times better than the Tevatron run II statistical 
expectations with 2 fb'^ (~ 0.03 on Fr WT\ and ~ 0.09 on Fq |58|). CDF and DO published 
first measurements of the W polarization in tt pairs based on Run I data [ !£, 20, 21 J. For exam- 
ple, CDF results are Fq = 0.91 ±0.37(stat)±0.13(syst) and Fr = 0.11 ±0.15. They are largely 
limited by statistical errors^. Preliminary studies started at Run II ll22ll . but statistical errors will 
remain large, even with the total integrated luminosity of 2 fb^^ ll59ll . 

The results of Table IT] were obtained assuming realistic uncertainties of 3% on the b-jet 
energy scale and 2 GeV on the top mass. More pessimistic assumptions (5% and 3 GeV) lead 
to an increase of the total systematic errors on Fi, Fq and Fr to 0.031, 0.018 and 0.013. On 
the contrary, more optimistic assumptions (1% and 1 GeV) lead to 0.018, 0.014 and 0.009. In 
all cases, the absolute error on Fr remains in the range 0.009-0.013 and that on Fq in the range 
0.015-0.018. The assumptions on systematic uncertainties have therefore a small impact, as- 
sessing the robustness of the results. 

All above results were obtained with a 10^-^cm^^s^^ luminosity. The luminosity may be 
2 ■ 10-^^cm^^s^^ at the LHC start. In this case, two scenarios are considered in ATLAS for 
the single electron trigger: increase of the single electron pj- cut at the trigger level from 20 
to 25 GeV (scenario 1), or even to 30 GeV (scenario 2). The complete study has been redone 
in the semileptonic channel for both scenarios, assuming the same hypothesis for each source 

few tens of both dileptonic and semileptonic events, with an integrated luminosity of 109 pb^^ (125 pb^') 
for CDF (DO). 
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of systematic uncertainty presented in Figure [131 The number of events will be multiplied by 
1.8 (1.6) for scenario 1 (scenario 2), while systematic errors remain almost unchanged. Con- 
sequently, the same precision will be achieved on the W polarization measurement. At high 
luminosity, lO^^cm^^s^^ a possible improvement can be to consider leptonic final states with 
J/\|/, in a similar way as what is performed for the top mass measurement iHTll . 

The forward-backward asymmetry, Afs, based on the angle between the charged lepton 
and the b-jet in the W rest frame, is often discussed in literature ll^ . It can be expressed in 
terms of Fi and Fr ifTTll : 

AFB = liFL-Fj,) (16) 

Taking the correlation between F^ and Fr into account, the following measurement on Afb can 
be extracted from the previous results: 

Afs = 0.227 ±0.003(stat)±0.016(syst) (17) 

Nevertheless, it does not provide any more information than the separate measurements of the 
ratios Fi and Fr. 

4.4 Sensitivity to new physics 

As already stated in the introduction, the search for anomalous (i.e. non Standard Model) inter- 
actions is one of the main motivations for top quark physics. The measurement of the W polar- 
ization provides a direct test of our understanding of the tWb vertex, responsible for practically 
all top quark decays in the Standard Model (SM). The deviations from the SM expectations 
induced by new physics contributions have been calculated in the framework of a few mod- 
els |l6Tl|62l|63l|64|- However, because of the great diversity of models beyond the SM (Super- 
symmetry, dynamical electroweak symmetry breaking models, extra dimensions, . . . ), it is also 
useful to study these possible new interactions in a model independent approach ITtI I^I^IFtI . 
The unknown dynamics can be parametrized with couplings representing the strength of effec- 
tive interactions, through the following Lagrangian [6J: 

L = ^W^bf{fi:PL + ffPR)t--^d^W^-ba^''{f^PL + f^PR)^ (18) 

where Pr/i = ^{1 ±y5), o-"^ = j['f,'f], 8 is the electroweak coupling constant and A is the 
energy scale to which the new physics becomes apparent (in the following, A =Mw is set to 
keep the notation used in the literature), /f and /f are vector- like couplings, whereas /2 and 
f2 are tensor-like couplings. This is the most general CP-conserving Lagrangian keeping only 
the leading (mass dimension 4, first term) and the next-to-leading (mass dimension 5, second 
term) effective operators in the low energy expansion. In the SM, the values of the couplings at 
tree level are /f = y,fc = 1 , /f = = /f = 0. 

It will first be shown how the measurement of the W polarization in top decay can probe 
these anomalous couplings (section l4.4.1l) and then a review of their existing direct and indirect 
experimental limits will be given ("section l4.4.2l) . 
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4.4.1 Probe of tWb anomalous couplings 



The contributions of ff', ff, and /j^ anomalous couplings to each fraction of helicity state, 
Fl, Fq and Fr, have been calculated at LO [ 6", "681 and at NLO [60|, NLO effects being small. 
A deviation of from 1 has not been considered in the following, as the W helicity is not 
sensitive to it. This can be easily understood, as is proportional to Vtb, which can not be 
directly measured with tt pairs but only with single top quarks. 

In the following, an independent deviation of each anomalous coupling, , /^^ and 
is assumed. Figure [T31 shows the variation of F^, Fq and Fr with these couplings. F^ and Fr 
depend quadratically on ff, whereas Fq remains unchanged (full lines). Similarly, the three 
fractions of helicity states are sensitive to in a quadratic way (dashed lines). In these two 
cases, the sign of the coupling can not be determined and the sensitivity will be lowered by the 
quadratic behavior. The last case is the most interesting (dash-dotted lines): Fq and Ft depend 
almost linearly on with a slope=0.7 [69|, while Fr is unchanged. Thus the sign of /j^ can be 
determined: Fq > Fq^ and F^ < fI^ (Fq < Fq^ and F^ > Ff'^) signs the presence of negative 
(positive) anomalous coupling The sensitivity is also higher than for ff and . The pre- 
cision to which Fi, Fq and Fr can be measured (Table |7l) sets the sensitivity to each anomalous 
coupling. It is represented by grey bands in Figure [151 Fr is the most sensitive observable to 
probe and , whereas Fq is better for 




Figure 15: Fq, Fi and Fr dependence on the anomalous couplings ff (full lines), (dashed 
lines) and f2 (dash-dotted lines), /f = 1 is assumed. The expected 1<3 uncertainties on Fq, F^ 
and Fr measured in tt pairs after one LHC year at low luminosity (Tabled are indicated with 
grey bands. 
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Figure [161 (full lines) shows the overall sensitivity (statistics+systematics) to each anoma- 
lous coupling that can be expected from the W polarization measurement with ti pairs at LHC. 
Dashed lines represent the statistical sensitivity only^. The corresponding 2a limits (statis- 
tics+systematics) are given in Table[Hl The best sensitivity is obtained on due to the presence 
of the large linear dependence. It is of the order of the deviations expected by models like the 
Minimal Super symmetric Standard Model (MSSM) or the Topcolor assisted Technicolor model 
(TC2) |,68J. 




-0.4 -0.2 0.2 0.4 -0.2 -0.1 0.1 0.2 -0.1 -0.05 0.05 0.1 



f,R f^L f^R 

Figure 16: Sensitivity in a to the anomalous couplings ff, f2 and f2 extracted from the W po- 
larization measurement in tt pairs after one LHC year at low luminosity (10 fb^^). Full lines 
indicate the overall (statistics+systematics) sensitivity, while dashed lines represent the statis- 
tical sensitivity only. 



4.4.2 Comparison with existing limits 

The only existing direct limits on the tWb anomalous couplings can be placed from Tevatron W 
polarization measurements in tt pairs, which are limited by the low statistics. As an example, 
the run I result Fr < 0.18 at 95% C.L. [.19J translates in /f < 0.8. From run II expectations, 
statistical sensitivities to /f ~ 0.5 and /^'^ ~ 0.3 ll68l at 95% C.L can be achieved. The single 
top, which has not been experimentally observed so far, can provide further constraints on the 
tWb anomalous couplings from its production rate and kinematic distributions. At LHC, the 
expected 2o limits are -0.052 < < 0.097 and -0.12 < /f < 0.13 COI assuming a 5% sys- 
tematic uncertainty, and a statistical sensitivity to /f ~ 0.06 f7T| with 100 fb^^, one LHC year 
at high luminosity. However, all these studies do not include any detector effect and detailed 
evaluation of systematic uncertainties. The related limits are summarized in Table [H for each 
coupling. 



recent study at NLO on forward-backward asymmeti-y Afb in tt pairs f60l indicates a 3a statistical sensi- 
tivity with 100 fb^' on ~ 0.06, /2 ^ 0.03 and f2 ^ 0.003. This is in good agreement with our 3a statistical 
sensitivity on /f ~ 0.17, /f ~ 0.08 and /f ~ 0.012 obtained with 10 fb^^ only. 
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Indirect limits on the tWb anomalous couplings have already been derived from precision 
measurements. The b ^ sy and b sl^l^ decays proceed via an electroweak radiative penguin 
process [72|. As they include a tWb vertex, an anomalous coupling will result in a change of the 
branching ratios. The related limits on anomalous couplings are stringent: as an example, /f 
has to be less than 0.004 ll73l at 95% C.L. As the tWb coupling appears also in loop in Z decays, 
electroweak measurements from LEP/SLC give other indirect limits, mainly competitive on 
All these limits are presented for each coupling in the last two lines of Table [SI However, they 
are indirect, SM-dependent, and scenarios can be envisaged where other contributions lead to 
cancellations that invalidate these bounds. 

To conclude, it is worth to notice that our expected sensitivity to the right-handed tensor- 
like coupling is a factor 2-3 better than the best limit. In any case, the W polarization 
measurement in tt pairs and the single top studies at LHC will be complementary to determine 
the structure of the tWb vertex as precisely as possible. 











LHC(10fb-^) 


0.30 


0.13 


0.04 


(Stat.-i- Syst.) 








?f, Tevatron (2 fb"^) 


0.5 


0.3 


0.3 


(Stat, only) 








single top, LHC (100 fb"^) 


0.06 


0.07 


0.13 


(Stat.-i- 5% Syst.) 








b sy.sl^l^, B-factories 


0.004 


0.005 


0.4 


(indirect) 








Z decay, LEP 






0.1 


(indirect) 









Table 8: 2a limits on anomalous couplings f f, and At LO, in the SM, these couplings are 
equal to zero. The first line presents our results extracted from the W polarization measurement 
in tt pairs after one LHC year at low luminosity (10 jb^^). Expected limits at the Tevatron and 
with single top at LHC are shown in the next two lines. Current indirect limits from B-factories 
and LEP data are presented in the last two lines. 
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5 Sensitivity to top quark polarization in tt events 



As demonstrated in the previous section, the W polarization measurement provides a direct 
probe of the top decay mechanism. Using the same events, it is also possible to test the tt 
production by measuring the top spin asymmetries, A and Aj). Similarly as for W polarization, 
we will explain the method used to extract these asymmetries (section lSTTT) . present a complete 
study of systematic uncertainties in both semileptonic and dileptonic tt channels (section ls!2b . 
give the results combining both channels (section l53t . and finally discuss the related sensitivity 
to physics beyond the Standard Model (section |5]4|). 



5.1 Measurement method 



Similarly to the W polarization analysis (section|4|), selection cuts distort the parton level angu- 
lar distributions. Therefore, expressions given in Equation (fT2b are no longer unbiased estima- 
tors of the spin correlation observables. To correct for this bias, a weight is applied on an event 
by event basis, allowing to get back, as much as possible, the original asymmetry. One weight 
is applied per spin correlation observable (A and Ajy) and per channel (semileptonic and dilep- 
tonic). The weighting functions are computed by fitting the selection efficiency in cosGi cos 62 
(cos4>) bins for A (Aj)). This is shown in Figure[l71 with a mean weight set to 1, which is the 
equivalent of Figure[51in the W polarization analysis. The four ratios are fitted by a polynomial 
function to extract smooth corrections. The choice of polynomial order has been tuned to get 
the best fit quality. The A corrections have a wider range with respect to the Ad ones, proving 
that A/) is less affected by selection cuts. The correction functions, computed on an independent 
data sample, are then applied event by event on the analysis samples. 



^ 2 





P1 1.288 
P2 -0.5220 
P3 -1.247 


-d 




pi 1.1 104 

P2 -0.3033 
P3 -0.3290 
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1 
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P3 



1 .0440 
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Figure 17: Ratio between the distributions o/cos0i cos 02 (top) and cos4> (bottom) after selec- 
tion cuts and at parton level for semileptonic (left) and dileptonic ( right) tt events. The full lines 
are the results of polynomial fits. 
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The correction functions are extracted with a Standard Model scenario. In case of de- 
viation from the SM, the kinematic distributions can be affected, and the correction functions 
will be changed. This is illustrated in Figure [TS] (left plots) for different A (top) and Ao (bot- 
tom) input values in the dileptonic channel^*^. For this purpose, different mixtures of events 
with/without spin correlation effects have been generated. Applying the SM correction func- 
tion to these samples will therefore not correct completely for the bias induced by the selection 
cuts. Figure [TSI (right plots) shows with open circles the measured asymmetries A (top) and 
Ad (bottom) as a function of their input values after applying the SM correction function. The 
measurement is clearly biased. As for Figure [TT] in the W polarization measurement, we pro- 
ceed iteratively to overcome this problem. The SM correction function is first used. Then, in 
case of deviation from SM expectations a new correction function is calculated with this new 
asymmetry, and applied. The process converges after a few iterations, as seen in Figure fTSl 




cos * Input 



Figure 18: Left: Correction functions for different A (top) and A^ (bottom) input values in the 
dileptonic channel. The case A = 0.422 andA^ = —0.290 corresponds to the SM functions of 
Figure\^ Right: A (top) and Ao (bottom) after selection cuts and correction as a function of 
their input values. The open circles correspond to the measurement using the SM correction. 
Full triangles, stars, squares and diamonds correspond to the first, second, third and fourth 
iteration, respectively (see text). The dotted line is y = x. 



'''This channel has been chosen to illustrate the method, as the bias due to the event reconstruction is more 
pronounced than in the semileptonic channel. 
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5.2 Systematic uncertainties 

The same sources of systematic uncertainties as for the W polarization study ("section l4!2b are 
considered: five related to the generation (Q-scale, structure function, ISR-FSR, Z^-fragmentation 
and hadronization scheme), three to the reconstruction (Z?-tagging, b-jet miscalibration, input top 
mass), the background normalization and the pile-up influence. In the semileptonic channel, the 
light jet miscalibration is also taken into account as the least energetic jet in the top rest frame is 
used as spin analyzer. A particular attention was paid to the proportion of gg and qq processes 
involved in the tt pair production, which directly impacts the spin correlation (see Figure |4l). 
To study separately this effect, samples with different proportions of gg/qq have been generated 
from 82%/18% to 90%/10% by steps of 2%. Small and smooth dependences are observed with 
a slope of 0.006 (0.004) per % of gg/qq variation for A (Ao). 

The results obtained on A and Ad for different Z^-tagging efficiencies, b-jet miscalibra- 
tions, top masses and pile-up levels are detailed in Figure [T^l Linear behaviors are observed 
in both channels. All systematic uncertainties are listed in Tableland illustrated in Figure EDI 
Generation and reconstruction sources contribute roughly in the same proportion to the total 
error. The dominant generation contributions come from the Q-scale, the structure function and 
the Z?-fragmentation, while the reconstruction systematics are dominated by the b-jet miscali- 
bration and the top mass uncertainty. The total systematic error for A is 2.5 times higher than 
for A/). This is because the angles are computed in the ti rest frame, more difficult to reconstruct 
than the top and anti-top rest frames separately. The lower systematics in the dileptonic channel 
are explained by the choice of two ideal spin analyzers (charged leptons). 



Source of uncertainty 


Semileptonic channel 


Dileptonic channel 




A 


Ad 


A 


Ad 


Generation 










2-scale 


0.029 


0.006 


0.011 


0.003 


Structure function 


0.033 


0.012 


0.008 


0.005 


ISR 


0.002 


0.001 


0.001 


0.001 


FSR 


0.023 


0.016 


0.005 


0.000 


Z^-fragmentation 


0.031 


0.018 


0.007 


0.004 


Hadronization scheme 


0.006 


0.008 


0.005 


0.003 


Reconstruction 










b-tSLgging (5%) 


0.016 


0.011 


0.001 


0.001 


b-jet miscalibration (3%) 


0.045 


0.012 


0.013 


0.003 


light-jet miscalibration (1%) 


0.009 


0.000 


0.000 


0.000 


Input top mass (2 GeV) 


0.028 


0.013 


0.009 


0.001 


Others 










S/B scale (10%) 


0.000 


0.000 


0.003 


0.004 


Pile-up (2.3 events) 


0.001 


0.005 


0.004 


0.003 


TOTAL 


0.081 


0.036 


0.024 


0.010 



Table 9: Summary of systematics on A and Ad in the semileptonic and dileptonic tt channels. 
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Figure 19: Measured A {left) and {right) in the semileptonic {black circles) and dileptonic 
{open circles) tt channels as a function of different parameters, see text for more details. Linear 
fits are superimposed in each case, and the corresponding slope is indicated. 
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Figure 20: Systematic uncertainties on A andAo in the semileptonic and dileptonic tt channels. 
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5.3 Results 



TablefTOlpresents the expected Standard Model results for A and A/? after one LHC year at low 
luminosity (lO^^cm^^s^^ 10 fb^^). In the semileptonic channel, the sensitivity is driven by 
the systematic uncertainties, which largely dominates the statistical ones, while both errors are 
comparable in the dileptonic channel. Combining the results of both channel studies, assuming 
a pessimistic 100% correlation of systematic errors, lead to the results shown in Table flOl right- 
most column. They allow to observe and measure the Standard Model spin correlation with a 
4% precision. 





Semileptonic (±stat±syst) 


Dileptonic (±stat±syst) 


Semilep-i-Dilep 


A 
Ad 


0.422 ±0.020 ±0.081 
-0.288 ±0.012 ±0.036 


0.404 ± 0.020 ± 0.024 
-0.290 ±0.011 ±0.010 


0.406 ± 0.014 ± 0.023 
-0.290 ± 0.008 ± 0.010 



Table 10: Standard Model results for spin correlation observables after one LHC year of data 
taking (lO^^cm^^s^^, 10 fb^^) in semileptonic and dileptonic tt channels. A combination of 
both results is presented in the last column. 

This result can be compared with the 40% precision expected from Tevatron run II with 
2 fb^^ neglecting the systematics. Experimentally, the tt spin correlation has never been 
observed. The DO experiment sets a lower limit on A with 6 dilepton events from run I 
(110 pb^^) |74|. This limit, A > —0.25 at 68% confidence level, can not be compared to the 
LHC values because the dominant production process at Tevatron is qq tt, and the Standard 
Model prediction is A = 0.88. 

The results of Table flOl were obtained assuming realistic uncertainties of 3% on the Z?-jet 
energy scale and 2 GeV on the top mass. More pessimistic assumptions (5% and 3 GeV) lead to 
an increase of the total systematic errors on A and A/) to 0.030 and 0.01 1. On the contrary, more 
optimistic assumptions (1% and 1 GeV) lead to 0.015 and 0.009. In all cases. Ad remains in 
the range 4%-5%. The assumptions on systematic uncertainties have therefore a small impact, 
assessing the robustness of this result. 

All above results were obtained with a 10^^cm~^s^^ luminosity. As already discussed in 
section H31 the luminosity may be 2 ■ lO^^cm^^s^^ at the LHC start and modify the pt electron 
cut at the trigger level. The complete study has been redone in the semileptonic channel for the 
two foreseen scenarios, assuming the same hypothesis for each source of systematic uncertainty 
presented in Figure l20l The number of events will be multiplied by 1.9 (1.8) for scenario 1 
(scenario 2), while systematic errors remain unchanged. Consequently, the same sensitivity 
will be achieved for A and A/) measurements. 
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5.4 Sensitivity to new physics 

As already stated in the introduction, a tt spin correlation observation would check that the 
top quark decays indeed as a quasi-free quark, i.e. in particular before hadronization can take 
place which could dilute the spin information. A measurement of the expected Standard Model 
spin correlation would test the top properties, with a left-handed coupling and a 1/2 spin. On 
one hand, this would allow to set an upper limit on its lifetime, directly linked to Cabibbo- 
Kobayashi-Maskawa matrix elements ll75l . On the other hand, this would offer a unique oppor- 
tunity to study a "bare" quark, free from long distance effects of QCD such as hadronization 
and confinement. 

Therefore, a possible deviation of tt spin correlation from the SM prediction will be a hint 
of new physics. Its measurement can be used to probe the presence of new interactions. For 
example, gtt anomalous couplings, linked to chromoelectric [ 8 1 and chromomagnetic [9 1 dipole 
moments which naturally arise in dynamical electroweak symmetry breaking models such as 
technicolor or topcolor, can affect the resultant ti spin correlation [76J. This is also the case 
in the presence of either a new heavy resonance in the tt production, such as a spin neutral 
Higgs boson lfT2ll (gg ^ H ^ tt), or spin 2 Kaluza- Klein (KK) gravitons lfT3l . As an example, 
in theories with large extra dimensions lITTI . the j'-channel mediated by graviton KK modes 
gives rise to characteristic spin configurations and angular distributions for outgoing particles, 
which reflect the spin-2 nature of the intermediate KK gravitons. With the sensitivity quoted 
in the previous section, a 5o deviation from the SM tt spin correlation can be observed if the 
fundamental scale of the extra dimensional theory is below 1 .5 TeV. 

New interactions in the decay can also affect the tt spin correlation. As an example fT5\, 
if a sufficiently light charged Higgs boson exists, such as in supersymmetric models, the decay 
t H^b can compete with the SM decay mode t W^b. As the charged Higgs decay to elec- 
trons and muons is largely suppressed, the deviation on the V7 -polarization measurement can be 
small. Contrarily, for m//+ < 150 GeV and at small tan (3 (< 2), the decay in two jets is favored, 
affecting the spin correlation in the semileptonic channel. As a result, with m//+ ~80 GeV, a 
5o deviation from the SM tt spin correlation can be observed if the branching ratio for top into 
charged Higgs plus Z^-quark is larger than 25%. 

6 Conclusions 

Because of its high mass, close to the electroweak symmetry breaking scale, the top quark is 
an ideal place to search for physics beyond the Standard Model. W polarization in top decay 
and top spin observables reflect in detail the interactions involved in top quark production and 
decay. Moreover, they can directly be inferred from the angular distributions of their respective 
decay products. Therefore, they give a good opportunity for precise tests of these interactions 
and are sensitive probes of new physics. Their precise measurements will be possible at the start 
of the LHC data taking, thanks to the very large sample of top events that will be accumulated. 
They will be complementary to Vth and cross section measurements, as NLO QCD corrections 
and theoretical errors are much smaller, of the order of 1%. 
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The ATLAS capability to measure the W polarization components Fq, Fl and Fr (for 
longitudinal, left-handed and right-handed helicity fractions) and tt spin asymmetries {A, Au) 
has been studied in the complementary semileptonic and dileptonic tt channels. The results of 
both channel studies have been combined. Leading-order Monte Carlo generators were used as 
well as a fast simulation of the detector. The clean signature of semileptonic tt events, a high 
statistics (around 100000 signal events after selection and reconstruction in one year at low 
luminosity, 10 fb~^) and a high signal over background ratio (more than 10) are the attractive 
features of this channel. In the dileptonic channel, the event topology reconstruction is compli- 
cated by the presence of two neutrinos in the final state, but feasible, the correct solution being 
found in 65% of the events. Even if the statistics and the signal over background ratio are lower 
than in the semileptonic channel, it is an attractive channel for the top spin asymmetry measure- 
ment, because the two charged leptons of the final state are the most powerful top spin analyzers. 

In both channels, selection cuts bias the measurements. A weighting method was set up 
to correct for it, and its robustness assessed. The sensitivity of the measurements is driven 
by the systematic uncertainties, which already dominates the statistical ones after one year at 
low luminosity. The main contributions to the total uncertainty come from the scale, the 
hadronization scheme, the FSR knowledge, the &-jet energy scale and the top mass. After one 
LHC year, the Standard Model parameter Fq can be measured with a 2% accuracy and Fr with 
a 1% precision, comparable to the expected precision on the top mass. Using the same selected 
events, the Standard Model top spin asymmetry can be measured with a precision around 4% 
with 10 fb~^ These results are robust against other hypothesis for systematic uncertainties and 
trigger scenarios. 

The sensitivity to physics beyond the Standard Model can be deduced from the above 
results. This has been studied in a model independent approach on the decay side by introduc- 
ing three tWb anomalous couplings, /f , /2 and /j^, which parametrize new physics. The best 
sensitivity, a 2a limit of 0.04, is obtained on /j^, which is better than indirect limits and expec- 
tations from other measurements. Finally, the sensitivity of the top spin measurement to new 
interactions such as a top decay to charged Higgs boson or new s-channels (heavy resonance, 
gravitons) in tt production have also been assessed. 
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